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ABSTRACT 

We perform a weak lensing and photometric study of the z = 0.42 supercluster 
MS0302+17 using deep / and V band images taken with the UH8K CCD mosaic 
camera at the CFHT. We use archival ROSAT HRI data to estimate fluxes, gas masses 
and, in one case, the binding mass of the three major clusters. We then use our CCD 
data to determine the optical richness and luminosities of the clusters and to map 
out the spatial distribution of the early type galaxies in the supercluster and in other 
foreground and background structures. We measure the gravitational shear from a 
sample of ~ 30, 000 faint background galaxies in the range 22 <^ mj <^ 26 and find 
this correlates strongly with that predicted from the early type galaxies if they trace 
the mass with M/L B ~ 250/i. We make 2-dimensional reconstructions of the mass 
surface density. These recover all of the major concentrations of galaxies and indicate 
that most of the supercluster mass, like the early type galaxies, is concentrated in 
the three X-ray clusters, and we obtain mean mass-to-light ratios for the clusters 
of M/L B ~ 260/i. Cross -correlation of the measured mass surface density with 
that predicted from the early type galaxy distribution shows a strong peak at zero 
lag (significant at the ~ 9-sigma level), and that at separations ^ 200/?r 1 kpc the 
early galaxies trace the mass very accurately. This conclusion is supported by 
cross -correlation in Fourier space; we see little evidence for any variation of Mj L or 
'bias' with scale, and from the longest wavelength modes with A = 1.5 — 6/i _1 Mpc 
we find M/L ~ (280 ± 40)/i, quite similar to that obtained for the cluster centers. We 
discuss the implication of these results for the cosmological density parameter. 
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1. Introduction 



Weak lensing, in the sense of the statistical distortion of the shapes of faint background 
galaxies, has now been measured for quite a number of clusters: [Tyson, Wenk, & Valdes 1990| ; 
Bonnet et al."T993|; gaEEnan et al. 199% [Bonnet, Mellier, & Fort 1994[ D ahle, Maddox, & Lilje 
1994 ; [Mellier et al. 1994[; port & Mellier 1994|; |Smail & Dickinson 1995|; Tyson & Fischer 
1995; [Fort et al. 1996|; jSeitz et al. 1996|; [Bower & Smail 1997|; |Smail etal. 1997|; F ischer 



et al. 1997; pischer & Tyson 1997[ and provides a powerful and direct measurement of the total 



mass distribution in clusters. Following the pioneering attempt of [Tyson et al. (1984[ ) using 



phot ographic plates, several groups dBrainerd, Blandford, & Smail 1996|; D ell' Antonio & Tyson 
1996; [Griffiths et al. 1996| ; |Hudson et al. 1998P have reported measurements of the 'galaxy-galaxy 
lensing' effect due to dark halos around galaxies, and there have also been attempts to estimate 
the shear due to large-scale structure (|Valdes, Jarvis, & Tyson 1983[ Mould et al. 1994[ S chneider 



et al. 1998). Most of these studies have been made with fairly small CCD detectors; this severely 
limits the distance out to which one can probe the cluster mass distribution and also limits the 
precision of galaxy-galaxy lensing and large-scale shear studies. 

Here we present a weak lensing analysis of deep / and V photometry of the field containing 
the z ~ 0.42 supercluster MS0302+17 taken with a large 8192x8192 pixel CCD mosaic camera 
(the UH8K; |Luppino et al. 1995| ) mounted behind the prime focus wide field corrector on the 
CFHT. The field of view of this camera on this telescope measures 0°.5 on a side, and greatly 
increases the range of accessible scales, and should yield improved precision for large-scale and 
galaxy-galaxy shear measurements. Our goal here is to explore the total (dark plus luminous) mass 
distribution in this supercluster field and see how this relates to the distribution of super-cluster 
galaxies and the X-ray emitting gas. 

The target field, centered roughly on RA = 3 h 5 m 24 s .O, DEC = 17°18'0".0, (J2000) 
contains three prominent clusters in a supercluster at z ~ 0.42. The first of these clusters to 
be found, CL0303+1706, was detected optically by pressler & Gunn (1992| ). It has a redshift 
of z = 0.418, a velocity dispersion a v = 912 km/s, and was the the target of an Einstein IPC 
pointed observation. This observation revealed two sources within ~ 1' of the position given by 
Dressier and Gunn, one of which, it is now apparent, coincides with densest and most massive 
concentration of galaxies in the cluster, and which has an X-ray luminosity L x — 3.0/i5 2 x 10 44 
erg/s. Analysis of the IPC image and optical follow up by the EMSS team ( pioia et al. 1990| ; 
Stocke et al. 1991| ; |Gioia & Luppino 1994fj revealed the presence of the two neighboring clusters 
MS0302+1659 and MS0302+1717 with redshifts z = 0.426,0.425 and X-ray luminosities 
Lx — 5.0 x 10 44 /i5Q 2 erg/s 4.3 x 10 44 /i^ 2 erg/s. The fluxes here are taken from the analysis of the 



IPC images by pabricant, Bautz, & McClintock (1994] ), who also obtained redshifts for a large 
number of galaxies in the the two X-ray brighter clusters, and obtained velocity dispersions of 
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and 921 and 821 km/s respectively. |Carlberg et al. (199"6j ) have reported a velocity dispersion of 



646 km/s for MS0302+1659 based on 27 redshifts and piling son et al. (1997[ ) have presented an 
extended catalog of redshifts for this cluster. Deep optical CCD images of MS0302+1659 were 
taken by [Mafhez et al. (199"2|) which revealed the presence of a pair of giant arcs and additional 



photometry was taken by |Giraud (1992| ) who found evidence for variability of features in one of 
the arcs. 

The angular separations of the clusters CL0303+1706, MS0302+1659 and MS0302+1717 
(which we will refer to as CL-E, CL-S, CL-N respectively) are 9 NE = 17'.64, 9 NS = 18'.49, 
and 9 SE = 14'. 04. Assuming an Einstein de Sitter universe, the angular diameter distance 
is D = 2cH 1 (l - (1 + z)- l / 2 )/(l + z) so D(z = 0.42) ~ 679/i" 1 Mpc and the physical 
scale is 3.29/iT 1 kpc/arcsec or 197/i 1 kpc/arcmin. For pure Hubble flow the line of sight 
physical separation is Ar" = (1 + z)~ 5 / 2 cAz/H or about 1.24A2/i _1 Gpc. Fabricant et 
al. found Az NS = 0.0009 ± 0.0012 from which we infer that the separation of the EMSS 
clusters in redshift space is ArNs — 3.8/i _1 Mpc, and is primarily transverse to the line 
of sight (Arj^g ~ 1.1 ± 1.5/i _1 Mpc, Ar^s — 3.6/i _1 Mpc). The Dressler-Gunn cluster 
E lies about Ar|j E ~ 9/i _1 Mpc in front of the two EMSS clusters. The parallel distance 
components may be somewhat distorted by peculiar motions. For an empty open cosmology 
D = 0.5cH \l - (1 + z)~ 2 ) so D(z = 0.42) ~ 756/i~ 1 Mpc and Ar'l = (1 + z)~ 2 Azc/H 
or about 1.49Az/i _1 Gpc, in which case the transverse and line of sight physical dimensions are 
larger than the EdS values by about 11% and 20% respectively. For an extreme A dominated 
model with Q A = 1, O m = 0, D = czH Q l /(l + z) and so D(z = 0.42) ~ 887/i" 1 Mpc and 
Ar" = (1 + z)~ l Azc/ Hq or about 2.1lAz/i _1 Gpc, in which case the transverse and line of sight 
physical dimensions are larger than the EdS values by about 30% and 69% respectively. The 
parallel distance components may be somewhat distorted by peculiar motions. 

This field was chosen for weak lensing as the high velocity dispersions and high X-ray 
luminosities of these clusters (along with the presence of arcs in one case) lead one to suspect that 
they might be potent lenses. Also, the whole system fits conveniently within the approximately 
30' x 30' field of view of the UH8K camera. This system provides an example of a quasi-linear 
system in the process of formation and it is of interest to see if weak lensing can detect, or 
constrain, the presence of ' bridges' of dark matter that can be expected between clusters (B ond, 
Kofman, & Pogosyan 1996). However, our main goal is to explore the correlation between 
mass and optical luminosity density to obtain estimates of the mass to luminosity ratio M/L, or 
equiv alently the inverse of the 'bias factor', (using the methodology of |Kaiser (199T|), S chneider 



(19985T and to see if this varies with scale. A source of systematic uncertainty in weak lensing 
mass estimates is the still largely unknown redshift distribution of the faint background galaxies. 
As a check we compare our mass estimates for the centers of the clusters with mass estimates 
from the velocity dispersions and, in one case, an estimate of the X-ray gas temperature. 
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The outline of the paper is as follows: In §^ we review the existing X-ray and optical 
spectroscopy results. We obtain improved fluxes and positions from archival ROSAT HRI 
observations. We make crude estimates of the central gas masses and, for CL-S estimate the 
central total mass from an ASCA temperature obtained by |Henry (1998| ) and compare with the 



optical velocity dispersion. In §|3] we explore the optical structure of the supercluster. We briefly 
describe the image processing; present a color image of the field; and show the size-magnitude 
distribution for the ~ 44, 000 objects detectable at > 6-sigma significance. We describe the 
somewhat non-standard photometric system we are using, and compute expected colors as a 
function of redshift and galaxy type. We plot the position, color and luminosities of the brighter 
galaxies, and we compute counts of galaxies as a function of magnitudes for regions around the 
clusters and for the complementary field region. In the absence of extensive redshift information 
it is hard to accurately measure the total optical luminosity distribution due to the confusing effect 
of foreground structures. If we restrict attention to the narrow band of V — I color occupied by 
the early type galaxies at the supercluster redshift this foreground clutter is effectively removed 
and we obtain a cleaner picture of the distribution of these galaxies and obtain estimates of the 
optical richness and luminosity functions for the three X-ray bright clusters. 



In §Q we present the weak lensing analysis. In §[44] we describe how we select the faint 
galaxies used to measure the shear and how we calibrate the relation between the observed image 
polarization and the gravitational shear. The calibration and also the correction for point spread 



function anisotropy is described more fully elsewhere ([Kaiser 1998|). Our shear analysis features 



an optimized weighting scheme, and we show how the weight is distributed in magnitude and on 
the size-magnitude plane. In §£T2] we predict the dimensionless surface density n{9) = S m /S crit 
assuming that early type galaxy light traces the mass with some a constant M/L. We do this both 
for a narrow slice around the supercluster redshift and also for a broader range of V — I color, 
corresponding to a range of redshifts 0.2 <; z ^ 0.6, and which includes a significant contribution 
from a foreground structure at z ~ 0.3. In §fO we correlate the measured shear with that predicted 



from the early type galaxy distribution. In § ftT4| we perform 2-dimensional mass reconstructions 
and compare these with the constant Mj L prediction. In § |4"3| we estimate central mass-to-light 



ratios for the three X-ray bright clusters, and we compare these with velocity dispersions and 



X-ray temperature information. In §[4^6] we compute the real space light-mass cross -correlation 
function to obtain a more precise estimate of the MJ L parameter and we compare the profile 
of the mass-light cross-correlation function and that of the luminosity auto-correlation function. 
We also perform the cross-correlation in Fourier space in order to see if there is any evidence for 
'scale dependent biasing' (i.e. variation of the M/L parameter with scale). 
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2. X-Ray Properties 

In this section we review the published X-ray results and obtain improved fluxes and 
positions from archival ROSAT HRI observations. We make crude estimates of the central gas 
masses and, for CL-S estimate the central total mass from an ASCA temperature obtained by 
Henry and compare with the optical velocity dispersion. 

The Einstein IPC image from the Einstein database is shown in the left hand panel of figure 
|T], with the circles indicating the locations of the clusters. Only the region where we have optical 
photometry coverage is shown. Analysis of these data by |Fabricant, Bautz, & McClintock (1993]) 



gave estimates of the luminosities for these clusters quoted in the Introduction. The identification 
of the target cluster was originally somewhat confused as the cluster is quite extended, and the 
location provided by Dressier and Gunn was displaced somewhat to the West from what is now 
seen to be the cluster center, and lay right between the two sources on the left with 5 ~ 17°19'. 
There is some indication of extended emission around CL-N, but the situation is somewhat 
confused by the shadowing of the support structure. 

There have been two further observations of the field with the ROSAT HRI which are now 
in the public domain. One of these was of ~ 3.5 x 10 4 s integration time and was centered on 
CL-S, but also contains CL-E, while the other of ~ 9 x 10 4 s integration contains all three clusters. 
These have higher resolution and higher astrometric precision than the IPC data. To analyze these 
data we used Steve Snowden's software ( {Snowden 1998| ) to generate a model for the particle 



background for each pointing, which we subtract, and also to generate a detailed exposure map 
which we divide into the subtracted image to correct for vignetting. We then mask out objects 
detected at > 4-sigma significance level and take the mean of what is left to obtain an estimate 
of the sky background, which we also subtract. The neutral hydrogen column density along the 
line of sight is Nh = 1 x 10 21 cm~ 2 ( Pickey & Lockman 1990D from which we find (assuming a 
Raymond-Smith spectrum with observed temperature of 3.43keV and 0.25 times solar metallicity) 
that one HRI count per second corresponds to an unabsorbed bolometric flux of 1.06 x 10~ 10 
erg/sec/cm 2 (for the often quoted passbands of 0.2 — 4.5keV and 0.5 — 2.0keV the corresponding 
conversion factors are 7.75 x 10~ n erg/sec/cm 2 and 3.79 x 10 -11 erg/sec/cm 2 ). A smoothed 
surface brightness image is shown in figure [2[ The unsmoothed HRI image shows that most of the 
other sources are still unresolved (and are therefore most probably not emission from cluster gas). 
The smoothed HRI image also shows the extended emission around CL-N, though now without 
the masking effect of the IPC ribs. CL-E is seen to be extended towards the NW, and at higher 
resolution appears to have a second core. This bimodality is also seem in the highly non-circular 
distribution of cluster galaxies so it would appear that this is most probably a merger event. The 
azimuthally averaged surface brightness profiles for CL-N, CL-S are shown in figure 0. CL-N 
appears to have a slope with index 7 ~ — 1 (where 1(6) oc 6 1 ) with significant emission extending 
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to 6 ~ 100 - 200". CL-S has a slightly steeper slope 7 ~ -1.4 within 9 ~ 50", outside of which 
radius the emission drops sharply and there is no statistically significant flux detectable at larger 
radii. This cluster is remarkably compact as compared to more typical bright X-ray clusters which 
have cores with essentially flat emission profiles extending to about this radius and with most of 
the gas residing at larger radii. 
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Fig. 1. — left hand panel shows the IPC X-ray image from the Einstein archive. The units of 
intensity are in IPC counts per second per square arcmin. The circles show the optical cluster 
centroids. Right hand panel shows a composite of the two HRI pointings, smoothed with a 35" 
Gaussian. The intensity units are HRI counts per second per (5") 2 pixel. 

Assuming spherical symmetry, the rest-frame surface brightness at impact parameter r is 

Irest(r) = ^Jdz e ff (v / ^T^) (1) 
where the Bremsstrahlung emissivity is eg = as^T 1 / 2 where as — 3 x 10~ 27 erg cm 3 s _1 Kr 1 / 2 



flSarazin 1988D - Approximating the integral as a sum over cubical 'voxels', the observed surface 
brightness for a pixel at radius i pixels from the cluster center in an image with pixel size AO is 

a uiA6 
4tt(1 + zf z - 



(2) 



with F{i) = eff(aiLUiiA9). We can invert the projection sum numerically to obtain 

F(i) = 4vr(l + zf(a Q uJiMy 1 I(i) where = EAV* 2 + J 2 )> and hence the gas 
density: 



n p (r) 



k ao^AOTWas 



1/2 



(3) 
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Fig. 2. — Upper panel shows the azimuthally averaged surface brightness profiles for the clusters 
CL-S (squares) and CL-N (circles). CL-E is not shown as it is clearly non-circular. Lower panel 
shows the gas density profile obtained by direct deprojection as described in the text. 
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The gas density profiles are shown in the lower panel of figure ^ from which we estimate gas masses 
within radius corresponding to 1'.5 of M s ~ 2.5 x lO 12 /i~ 5/2 M , M N ~ 2.8 x lO 12 /^ 5/2 M . 
I H 

Henry (1998) has obtained an ASCA spectrum for CL-S from which he finds a rest- 
frame temperature of 4.6 ± 0.8keV. Assuming isothermality, the total mass interior to r is 
M(r) = —(kTr /Gfim p )dlnn/dlnr, and using a logarithmic slope of d In n/d In r = —1.5 and 
a mean molecular weight fi = 0.6 we find M(< r = 167/?r 1 kpc) = 4.25 x \Q Vi h~ l M Q , this 
radius corresponding to the angular scale of 50". Combining these results, we find a gas fraction 
for CL-S of M gas /M to t — 0.06/i~ 3 / 2 . Both the gas and total mass estimates should be considered 
to be lower bounds, as there may be emission at larger radii which we cannot detect and one 
could also imagine that the cluster has an extended dark matter halo which these central X-ray 
observations cannot constrain. It is interesting to compare (kT/Gfirrip) 1 / 2 ~ 854km/s with the 
observed velocity dispersion of 921km/s. As we shall see, the galaxies in the cluster are, like 
the observable gas, highly concentrated and have a very similar scale length, so it is reassuring 
that the kinetic energy per unit mass is essentially the same for the gas and for the galaxies. The 
high temperature and very small size of the cluster indicate a very short sound crossing time 
t = r/c s ~ 2 x 10 8 yr. 

The inferred gas density in the very center of CL-S (at around 15" or about 50/i _1 kpc say) 
is n ~ 1.4 x 10~ 2 /i _1 / 2 cm^ 3 , implying a cooling rate of about T ~ 1.4 x 10~ 18 /i~ 1 / 2 s _1 or a 
cooling time t coo i ~ 2.2 x 10 9 /i 1 / 2 yr, as compared to the age of the Universe at z = 0.42 of 
t ~ 3.9 x 10 9 /i _1 yr for O = 1, or t ~ 7.1 x 10 9 /i _1 y r for an open £1 ~ cosmology, so the ratio 
of the cooling time to the age of the Universe t coo \/t(z = 0.42) is in the range (0.3 — 0.5)h 3 ^ 2 
so we are dealing with a strong, but still quasi-static, cooling flow. We estimate the gas loss rate 
as M ~ 100 — 2OOM /yr. At a radius of 50" the density is an order of magnitude lower, so the 
cooling time is an order of magnitude larger and, for h = 0.7, t coo i/t(z = 0.42) lies in the range 
1 .7-3.0 so the cooling is somewhat larger than the age of the Universe. However, for low Q, and in 
the absence of heat sources, much of the gas we see will by now have dropped out of the cooling 
flow. 

The HRI derived positions, count rates, fluxes luminosities of the three clusters are 
summarized in ta ble jl]. The agreement between our luminosities and those of Fabricant, Bautz, & 
McClintock (1994) is generally very good. For CL-S, for instance, and for the 0.2 — 4.5keV band, 
we find L = 4.9 ± 0.4 x /^ 2 10 44 erg/s whereas they find L = 5.0 ± 0.8 x /i^lO 44 erg/s. 



3. Optical Photometry 
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Cluster 


RA(J2000) 


DEC(J2000) 


r [cts/s] 


/boi [erg/s/cm 2 ] 


hl L ho i [erg/s] 


CL-S 


3 ft 05 m 31 s .49 


17°10'16".32 


7.1 ±0.5 x 10" 3 


7.6±0.6" 13 


6.8 ± 0.5 x 10 44 


CL-E 


3 h 06 m 18 s .98 


17°18'33".91 


3.8 ±0.7 x lO- 3 


4.0±0.8~ 13 


3.6 ±0.7 x 10 44 


CL-N 


3 h 05 m 17 s .81 


17°28'37".64 


5.6 ±0.6 x 10" 3 


6.0±0.7" 13 


5.4 ±0.6 x 10 44 



Table 1: Cluster X-ray Properties. Positions are locations of peaks of the smoothed HRI image. 
Count rates, fluxes and luminosities are for 1'.5 radius apertures centered on the given position, 
except for CL-E where we include the flux in a second aperture centered on the secondary feature 
to the NW of the main cluster. 



3.1. Data Acquisition and Processing. 

In September 1995, we obtained deep images of the MS0302 supercluster using the UH8K 
CCD camera mounted at the prime focus of the CFHT. At an image scale of ~ 0".207/pixel, the 
8192x8192 pixel CCD mosaic spans a field of view of ~29' on a side. The images were taken in 
superb observing conditions: photometric sky with ~ 0".6 seeing. Multiple 15 minute exposures 
were taken, with the telescope shifted slightly (~ 20 — 30") between exposures, in order to 
facilitate fiat fielding and the construction of a seamless final image. In all, 1 1 I-band images 
and 6 V-band images were used to give final exposure times of the summed images of 9900s and 
5400s respectively. 

The i mage processing was rather involved, and is described more fully elsewhere (K aiser 
et al. 1988) but essentially consists of four phases; pre-processing of original images; registration 
to obtain the transformation from chip to celestial coordinates; warping and averaging to produce 
the final images (and also auxiliary files describing the variation of photon counting noise, point 
spread function shape etc); and finally the generation of catalogs of objects. 

In the pre-processing we first generated a median sky image or 'super-fiat' for each chip and 
each passband and divided this into the data images. The super- fiat was also used to identify bad 
columns, traps, and other defects, which were flagged as unreliable. To further suppress residual 
variations in the sky background level we then subtracted from the images a smooth local sky 
background estimate formed from the depths of the minima of a slightly smoothed version of the 
super-flat divided data images. 

In the registration phase we took the locations of stars measured on the data chips, along 
with stars with known accurate celestial coordinates from the USNOA catalog ( pVIonet 1998| ) and 



solved in a least square manner for a set of low order polynomials — one for each of the 2K x 4K 
images — which describe the mapping from rectilinear chip coordinates to celestial coordinates. 
We computed a sequence of progressive refinements to the transformation. The final solution 
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gave transformations with a precision of about 5 milli-arcsec, which accurately correct for field 
distortion, atmospheric refraction, and the layout of the chips in the detector, and also any other 
distortions produced by inhomogeneities of the filters and/or mechanical distortion of the detector. 

In the image mapping phase we applied the spatial transformations to generate a quilt of 
slightly overlapping images covering a planar projection of the celestial sphere. We chose the 
'orthographic', shape-preserving projection ( preisen & Calabretta 1995] ) with tangent point at 



(a, <5)j2ooo = (3 h 05 m 24 s .O, 17°18'0".0) and rotation angle of 2°.635 (this being chosen to render 
the star trails aligned with the vertical axis of the final image to facilitate masking of these 
artifacts). Cosmic ray pixels and satellite trails etc. were identified by comparing each image 
section with the median of a stack of images taken in the same passband, and were flagged as 
unreliable. The final images were made by simple averaging of the non-flagged pixels in the stack 
of warped images. The combination of the mosaic chip geometry layout, the pattern of shifts 
for the exposures, and the presence of unreliable data of various types resulted in a somewhat 
non-uniform sky noise level, which we kept track of by making an auxiliary image of the sky 
noise. The FWHM of stars in these stacked images were approximately 0".60. We also generated 
a detailed model for how the point spread function (psf) shape varies across the final image (this is 
quite complicated as the psf shape changes discontinuously across chip boundaries) and we apply 
this below to correct the measured galaxy shapes to what they would have been if measured with 
a telescope with perfectly circular psf. 

A color image generated from the summed V and I band exposures is shown in figure |3|. 
The Dressler-Gunn cluster CL-E is clearly visible close to the Eastern edge of the field as an 
agglomeration of early-type galaxies with very similar color. The two X-ray selected clusters 
can be seen fairly clearly, and the giant arcs in the Southern cluster are apparent. Both of these 
clusters contain galaxies with the same color as in CL-E, but in both cases there are also bright 
bluer galaxies in the vicinity, suggesting that there is some foreground contamination. 

After registering and summing the images we ran our 'hierarchical peak finding algorithm' 
to detect objects ( |Kaiser, Squires, & Broadhurst 1995| , |Kaiser et al. 1988| ) on each of the I- and 



V-band images and then ran our aperture photometry and shapes analysis routines to compute 
magnitudes, radii and other properties of the objects. For each of the I, V catalogs we also 
estimated magnitudes using the summed image for the other passband (but with the same center 
and aperture radius) in order to provide accurate color information. We also constructed a 
combined T+V catalog in which, if an object was detected in both V and I then we include 
only the most significant detection, whereas if the object was only detected in one passband 
then the information for that filter only is included. This is the primary catalog we use in the 
analysis below. To remove spurious detections from the diffraction spikes around bright stars 
and trails caused by poor charge transfer for saturated pixels, we generated a mask consisting of 
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MS0302+1 7 




03 h 06 m 20 s 06 m 00 s 40 s 20 s 05 m 00 s 40 s 20 s 

RA [J2000] 

Fig. 3. — Color image generated from the summed V and I band images. This is a heavily reduced 
(de-magnified) image of the field derived from the final averaged images (which have a re-sampled 
pixel scale of 0".15 and cover 12, 000 2 pixels). In order to show both the faint background galaxies 
and the much brighter foreground galaxies, the brightness of each pixel here was taken to be 
proportional to the logarithm of the sum of the two V and I-band images, and the individual 
intensities are such that the red and blue signals are in the same proportion as in the linear I, V 
images (with a constant multiplicative scaling applied to the blue value) and the green signal is the 
mean of the red and blue signals. Thus this is not a true 3-color image (which would look much 
whiter), but is useful to distinguish the different galaxy types present. 
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the rectangles shown in figure [TT|, and objects lying within these rectangles were not used in the 
following analysis. The combined catalog contains about 44,000 objects detected at > 6-sigma 
significance. The total unmasked area is A = 0.165 square degrees, so the number density of 
these objects is about 74 per square arcmin. The size-magnitude distribution for the combined 
catalog is shown in figure 0, from which it can be seen that the 6-sigma significance threshold of 
the catalog corresponds to mj ~ 25.5 for point-like objects — though with incompleteness setting 
in at brighter magnitudes for extended objects. 

SIZE— MAGNITUDE DIAGRAM 
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Fig. 4. — Size-magnitude diagram for the combined sample. The vertical stellar locus is clearly 
seen, and stars can clearly be separated from galaxies for mj £ 23. At fainter magnitudes the stars 
and galaxies merge, but the stars are then a small contamination. 



3.2. Photometric System and Galactic Extinction 

The photometric system used here, whose passbands we denote by B', V, R', I', and which 
differ som ewhat from the standard B, V, R, I passbands, is described in detail elsewhere (K aiser 
et al. 1988). There we compute the redshift and SED dependent A;-terms in the transformation 
from UH8K apparent magnitude m F / to standard passband absolute magnitude M F : 

M F — cp = mpi — cp + kpF'(z; type) — 5 log 10 (D;/10pc) (4) 

Here cp is the magnitude of Vega in passband F. We also show that one can 
form synthetic magnitudes as a linear combination of our UH8K I', V magnitudes 
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m Fz = a Fz (m r - cj) + (1 - a Fz )(m v > - c v ) 
standard rest-frame absolute magnitudes are 



for F = B, V . . ., whose transformation to 



M F = m Fz + k Fz (z) - 5 log 10 (A/10pc) - c F 



(5) 



and where the coefficients a Fz have been chosen so that the k term is essentially independent 
of spectral type for galaxies at the redshift of the clusters z = 0.42. For an Einstein - de Sitter 
Universe, and with h = 1, the distance modulus is 5 log 10 (A/10pc) = 40.68. 

The apparent magnitudes here have been corrected for galactic extinction using 
the [Schlegel, Finkbeiner, & Davis (1998] ) IRAS-based prediction. Their maps give 
E(B - V) = 0.125,0.134,0.113 for the clusters N,E,S respectively. These are similar, 
and we adopt an average value of E(B — V) = 0.124, from which we obtain corrections of 
Am/ = 0.243, Am v = 0.402. 

In the lensing analysis we shall use the distribution of bright galaxies to predict the projected 
dimensionless surface density and thereby the gravitational shear. In doing this it is important to 
discriminate between structures at low and high redshift (as the lensing signal derives largely from 
the latter). A useful discriminant, especially for early type galaxies, is the V — V color, which is 
plotted as a function of redshift for various galaxy types in figure [5[ 

V - I COLOR vs REDSHIFT 




Fig. 5. — Galaxy color vs redshift for the UH8K system and using broad band SED's from 
Coleman, Wu, & Weedman (1980| ) . Symbols encode the galaxy type as follows: EO = stars; 
Scd = squares; Sbc = triangles; Im = circles. 



- 16- 



3.3. Bright Galaxy Properties and Cluster Luminosities 

We will now explore the properties of the brighter galaxies which define the supercluster 
and estimate the optical luminosities of the clusters. In the absence of redshift information this is 
rather difficult because of the presence of foreground structures. We can however obtain fairly 
accurate estimates of the individual cluster luminosities for the early-type cluster galaxies as the 
noise due to foreground and background contamination is then greatly reduced. 

In figure || we plot the locations, colors, and fluxes of the bright galaxies in the field. At 
the center of each circle drawn around the X-ray positions resides a dense concentration of very 
red galaxies with colors as expected for elliptical galaxies at z = 0.42. In addition to the red 
elliptical galaxies, a large number of bluer and brighter galaxies can be seen. One conspicuous 
concentration lies 2' — 3' to the East of CL-N, and, to judge from the colors and fluxes of the 
galaxies, lies at z ~ 0.2. Another, more diffuse concentration, with colors and luminosities 
indicative of early type galaxies at z ~ 0.3 can be seen extending Westward of CL-S. We have 
obtained redshifts for two of the bright galaxies in this complex and find that they do indeed lie at 
z ~ 0.3. 

We can obtain an unbiased, if somewhat noisy, estimate of the abundance of galaxies in 
general in the three X-ray clusters by computing the counts of all galaxies (the number of galaxies 
per solid angle per magnitude interval) within the cluster-centered circles and then subtracting the 
'background' counts obtained for the complementary region. The counts as a function of /-band 
apparent magnitude are shown in figure [7|. In computing the counts we allow for the somewhat 
irregular survey geometry caused by our masking of bright stars (the mask is shown below in 
figure [11]). To estimate the unmasked areas we count the number of objects from a densely 
sampled random Poissonian distribution of points which was masked in exactly the same way as 
the galaxy catalogs. The error bars in figure [TTjare poissonian only, and the real fluctuations are 
somewhat larger. The excess of galaxy counts due to the clusters is clearly visible. Similar results 
are obtained for the synthetic rest-frame 5-magnitude m Bz ~ (m,/ + m v )/2, though these start to 
become incomplete at somewhat brighter magnitudes. 

The color magnitude diagrams shown in figure |8| reveal a distinct excess of galaxies in a 
narrow band of color around V — I = 2.2 in the cluster regions, as predicted for early type galaxies 
at z = 0.42 from figure |5|. By selecting only galaxies within a narrow strip in color-magnitude 
space straddling this sequence we can obtain fairly precise estimates of the early-type luminosity 
functions for the three clusters, and for the clusters taken together, since the fluctuations in the 
background estimates are greatly reduced. The results are shown in figure ^ where we plot the 
spatial distribution of galaxies in a ~ 0.5 mag wide strip (corresponding to a range of redshifts 



Az ~ 0.1) and in figure|10l where we plot the excess counts of the color selected galaxies 



(corrected for field contamination) AiV = A c x (n c — n/), where A c is the area of the cluster 
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m, < 21 
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Fig. 6. — Plot of all galaxies brighter than I = 21 where the area of each circle is proportional 
to the flux density of the galaxy, and the shade of the disc encodes the V — I color, ranging from 
V — I — l for a white disc to V — I = 3 for a black disc. The legend in the upper right corner shows 
the image that would be produced by an galaxy (M B = —19.68) of various types at redshifts 
0.2, 0.3, 0.4. Stars have been removed. The three main galaxy concentrations are apparent, as is an 
extension to the West from southern EMSS cluster. The large circles have radius 3' (1'.5 for CL-S) 
and are centered on the X-ray cluster positions. 
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CLUSTER AND FIELD COUNTS 




m. 



Fig. 7. — Counts of galaxies for the 'cluster' regions defined by the circles in figure and for the 
complementary region. 

circles, as a function of M B = mj/ — cp + c B + k BI i(z = 0.42) — DM(z). Also shown is the 
luminosity function of form similar to those found by Efstathiou, Ellis, & Peterson (1988D and 



ILoveday et al. (1992Q : <&{L B ) = <^*(L/L,) a exp(-L/L*) with M* = -19.68 and a = -1.0 and 



with amplitude scaled to fit the counts. This standard, low-redshift, field galaxy luminosity 
function seems to provide a reasonable fit to the cluster count excess down to magnitudes well 
below the knee, and so integrating the luminosity functions should give a good estimate of the 
total red cluster light. 

To obtain estimates for the individual cluster luminosities we fit the excess counts to the form 
above to determine as 

$* = E AN f (Mb) Ax 2 / E f(M B ) 2 /a 2 (6) 

where f(M B ) = exp(— L/L*) and a is the rms fluctuation in AN due to Poisson 

counting statistics. We then integrate the luminosity functions to obtain the total excess 
luminosities within the circles defining the cluster subsample. These are tabulated in table ^. 
These luminosity functions are only for the color selected subsample. Performing the same 
exercise to compute the total luminosity function using our synthetic 5-band magnitudes we find 
noisier results, but conclude that about 70% of the total excess light within our cluster apertures 



comes from early type galaxies, in rough agreement with [Fabricant, Bautz, & McClintock (1994| ) 
who found that about 25% of the galaxies in CL-N for instance were blue with 'Balmer-line' 
spectra. 
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cluster 


0* 


LIU 


h 2 L/L Q 


E+N+S 


59.27 


128.70 


1.49e+12 


CL-S 


6.34 


13.77 


1.60e+ll 


CL-E 


30.30 


65.81 


7.63e+ll 


CL-N 


17.87 


38.81 


4.50e+ll 



Table 2: Table of cluster luminosities. 
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COLOR SELECTED 
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Fig. 9. — Early type galaxy sample. The spatial distribution of galaxies with colors 1.9 < V — I < 
2.4 and / < 23 chosen to select early type galaxies at around the redshift of the super-cluster. The 
area of each circle is proportional to the flux. The range of colors here corresponds to a depth 
Az ~ 0.1, so we would expect to find a substantial contamination by 'field' galaxies in addition to 
the super-cluster galaxies. 
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COMBINED CLUSTER LF 
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Fig. 10. — Rest frame £>-band luminosity functions for the three X-ray clusters and for the three 
clusters combined. The counts here have been corrected for field contamination and the lines are 
fits to a Schechter function with M*, a fixed at the values obtained for low redshift field samples. 
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The standard definition of cluster richness (|Abell 195%|) is the number of cluster galaxies in 



the range two magnitudes below the third brightest, with richness class R = 1 corresponding to 
50-80 galaxies etc. It is clear from inspection of figure [FT] that these clusters are not very rich by 
this measure. Both CL-N, CL-E are borderline richness class 1 and CL-S is a very poor cluster. 

The southern cluster CL-S is particularly compact; most of the elliptical galaxies lie within 
about 1' (or about 200/?r 1 kpc) of the X-ray location. This is almost an order of magnitude smaller 
than the conventional Abell radius of 1500/i _1 kpc. This is consistent with the compact X-ray 
appearance for this cluster, but also somewhat surprising. According to the standard hierarchical 
clustering picture, clusters which turn around and virialize at time t have a density ~ 170 times the 
critical density at that age: p ~ 170p crit (t) = 170/(67rGt 2 ). With r ~ i^Abeii and cr ~ lOOOkm/s, 
as typical of massive clusters, we obtain a collapse time t ~ 10 10 /i _1 yrs, similar to the present 
age of the universe and therefore in reasonable agreement with the theoretical expectation. With 
a = 920km/s and r ~ 200/i _1 kpc however, we obtain a collapse time t ~ 2 x 10 9 yr, which is 
much less than the age of the universe at z = 0.42. If this cluster is virialized and did indeed form 
by something like the spherical collapse model then it did so a very long time ago at z ~ 3 — 5 
and little visible matter has accreted onto it since then. 



4. Shear Analysis 

4.1. Background Galaxy Selection and Shear Measurement 

We now describe how we select the faint galaxies which we use for measuring the shear, how 
we correct for psf anisotropy and calibrate the effect of seeing. We also describe how we construct 
a weighting scheme to combine the signal from galaxies of a range of sizes and luminosities to 
give a reconstruction with optimal signal to noise. 

Stars are easily visible to rrij ~ 23 (figure ^) and we remove these and also the tail of faint 
objects with estimates half-light radii less than 0.3". Some close pairs of galaxies get detected as 
single objects, and these were filtered by applying a cut in ellipticity (see below) which we set 
at e = 0.7. To obtain a sample which is distant and minimally contaminated by cluster galaxies 
we make a cut at v = 100 to remove bright objects. The resulting catalog contains about 30, 000 
objects, with a surface density of about 50 per square arcmin. The spatial distribution of this faint 
galaxy sample is shown in figure [□], which confirms that they are very uniformly distributed 
across the field. 



The shear measurement is described in more detail elsewhere ( |Kaiser 1988| ). It is similar in 
spirit to the [Kaiser, Squires, & Broadhurst (1995| ) analysis in that we compute response functions 
which tell us how the shape statistics respond to psf anisotropy and to gravitational shear. The 
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FAINT GALAXY SAMPLE 




Fig. 11. — Spatial distribution of the faint galaxy sample. They are clearly distributed very 
uniformly across the sky as one would expect if they are very distant. The rectangles were placed 
to mask out diffraction spikes and trails from bright stars. 
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fundamental shape statistic we use is the flux normalized weighted second central moment of the 
measured surface brightness f(r): 

J d 2 rnr m w(r)f(r) 
J d 2 rf(r) 

where w(r) is a Gaussian profile weight function with scale matched to the size of the object. 
KSB computed the response of the ellipticity or 'polarization' defined as e a = M a i m qi m / 'q pp , 
where the pair of 2 x 2 matrices are M = {{1, 0}, {0, -1}}, M x = {{0, 1}, {1, 0}}, to both 
psf anisotrop ies and, in the limit of well resolved galaxies, to gravitational shear. Luppino & 
Kaiser (1997) generalized this to obtain the post-seeing shear response function P 7 such that 
(e Q ) = laPy and which serves to calibrate the relation between the observable polarization and 
the gravitational shear 7 Q = (l/2)M aJm $j m where $j TO is the projected gravitational potential. 

In |Luppino & Kaiser (1997D the psf anisotropy was modeled as a low order polynomial in 
position on the image. Such a model is not adequate for observations such as those considered 
here taken with a mosaic camera where one finds that the psf varies smoothly across each chip, but 
changes discontinuously as one passes across chip boundaries resulting in a complicated pattern 
of psf anisotropies on the final averaged image. Since we know which source images contribute 
to a given galaxy image, we can compute an average of the psf anisotropy response function 
measured from the individual images modeling these as a smooth polynomial functions over 
each of the source images. This would not be quite correct however, since the response function 
for the ellipticity, which is a non-linear function of the surface brightness, is itself non-linear, 
so the averaged response function does not, in general, correctly describe the response of an 
averaged galaxy image. What we do instead is to compute the response functions for the flux 
normalized q im moments themselves. Since psf anisotropy does not affect the net flux appearing 
in the denominator (we assume that our fluxes are effectively total fluxes), the psf response is 
a linear function of the sky brightness /(r) and can therefore be averaged to give the correct 
response for the q\ m which we use to correct the qi m values to what would have been measured by 
a telescope with a psf of same size and overall radial profile as the real psf but with no anisotropy. 
We then form a polarization e a as above from the anisotropy-corrected moments, and compute 
the post-seeing shear response function P 7 , such that, on average, the polarization e a induced by 
a gravitational shear 7 Q is 

(e«) = P 7 7 a (8) 

We do not attempt to calibrate the galaxies individually; this tends to be noisy and introduces 
non-linear effects. Instead we split the galaxies up into a set of discrete subsets by a coarse 
binning in the significance - half light radius plane (u, r h ), and compute an average P 7 for all the 
galaxies in each bin and thereby compute shear estimates %, = e a / (P 7 ). Having split the galaxies 
up into classes by size and magnitude in this way we can then compute an optimal weight as a 
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function of the bin. To do this we assume that all of the galaxies are at large, or at least similar, 
distance from the lens and therefore feel essentially the same shear, and also that the shear signal 
is much less that the random shear noise due to intrinsic shape and measurement errors combined. 
We compute the variance in shear estimates (7n 0ise ) — (e 2 ) / (-P 7 ) 2 for each bin and assign weights 
inversely proportional to the variance and such that the sum over galaxies of the weights is equal 
to the number of galaxies (so the result is unbiased). With this scheme, faint and/or very small 
galaxies which are relatively noisy get assigned lower weights. 

In figure [12| we show the distribution of weight in the size magnitude plane and we also show 
the distribution of weight as a function of magnitude. The rms shear (per component) in the final 
catalog is er 7 ~ 0.32 so we should be able to measure shear to a high precision. For instance, for 
measurement of the net shear we expect the precision to be a 7 / yJN^s ^ 0.002 

2-D WEIGHT DISTRIBUTION WEIGHT DISTRIBUTION 




Fig. 12. — Left panel shows the distribution of weight as a function of size and magnitude. The 
edges of the coarse bins can be seen as discontinuities in the weight density. Right panel shows the 
distribution as a function of magnitude. Roughly half of our weight comes from galaxies brighter 
than mi = 23.5, where the redshift distribution is fairly well established flCowie et al. 1996| ), but 
half comes from fainter magnitudes where there is only very poor sampling as yet. 

The largely unknown redshift distribution for galaxies at this range of magnitudes and sizes 
is a major source of systematic uncertainty in weak lensing mass estimation. For nij ^ 23 there 
are reasonably complete statistical samples flCowie et al. 1996| ) but at fainter magnitudes things 
become progressively uncertain. For sources at a single redshift z s , the dimensionless surface 
density k that we measure for a lens of physical mass surface density S m is k = S m /S crit (-2 S ) 
where 



J crit 



(9) 



4ttGA D ls 

where Di s is the angular diameter distance relating physical distance on the source plane to angles 
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at the lens and where the factor (3 = D ls / D s is the distortion strength relative to that of sources 
at infinite distance. With sources at a range of distances what we measure is the physical surface 
density times an average inverse critical density, or k = (E m /E^ it )(/3). The critical surface 
density and the (3 parameter are shown, for a lens redshift z\ = 0.42 in figure [TJ. 



INVERSE CRITICAL SURFACE DENSITY 



BETA PARAMETER 





Fig. 13. — Left panel shows the inverse critical surface density as a function of source redshift for 
a lens at Zi = 0.42. The physical surface density of a lens is obtained by dividing k by a weighted 
average inverse critical surface density taken over the distribution of background galaxies. The 
panel on the right shows the (3 parameter; for a lens of given velocity dispersion, X-ray temperature, 
the predicted shear is proportional to (3. The three lines are for EdS (solid), empty open (dot-dash) 
and A dominated cosmology (dotted). 



Using the photometric redshifts in the Hubble Deep Field from ( [Sawicki, Lin, & Yee 1997D 



and assuming an empty open universe gives ((3) ~ 0.46. Using an alternative, and more recent, 
HDF photo-z catalog calculated flFern'andez-Soto, Lanzetta, & Yahil 1998[ ) we find similar results. 
Both catalogs contain about 1000 galaxies. The {(3) value was calculated by transforming HDF 
-^F8i4W AB magnitudes to Cousins /-band magnitudes according to the WFPC2 Handbook and 
HDF web-page, and corresponds to that for a single screen of background galaxies at z s ~ 1.0. 

There is considerable uncertainty in this value. There are some questions about the reliability 
of photometric redshifts, the calculation above does not incorporate the weighting as a function of 
size, and the cosmological background model may be inappropriate. Another result which casts 
some doubt on the detailed accuracy of the distribution of photometric redshifts is that they are 



somewhat hard to reconcile with weak lensing results for z ~ 1 clusters ( [Luppino & Kaiser 1997| ; 



CJlowe et al. T39"%D- In what follows we shall adopt a value of S crit = 4.33 x lO 15 /iM /Mpc\ 



slightly higher than that derived from the photometric redshifts, and corresponding to a single 
screen redshift of z s ~ 1.5, 1.2; 0.75 for EdS, empty open and A dominated cosmologies 
respectively. 
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4.2. Predicted Dimensionless Surface Density 

We now compute the dimensionless surface density k that would be generated by the 
optical structures assuming that galaxies are good tracers of the total mass with some constant 
mass-to-light ratio We know of course that on mass scales of the order of tens of kpc the dark 
matter distribution is much more extended than the light, but on larger scales this 'what you see is 
what you get' picture remains an interesting and viable hypothesis. 

If all of the structures lay at the same lens-redshift zi this would be straightforward since we 
then have 

k=^- = (M/L)Ei/E crit (^) (10) 

^crit 

where £/ is the rest frame luminosity surface density, which we can compute from the observed 
surface brightness. 

This would be adequate to compute the contribution just from the z = 0.42 supercluster, but, 
as we have seen, there are other structures at different redshifts that we would like to include, and 
one cannot then use this simple procedure because it would cause one to grossly overestimate the 
effect of foreground structures. This is because, for a lens with given rest-frame surface luminosity 
and mass surface density, the observed surface brightness falls very rapidly with increasing lens 
redshift (roughly as (1 + z)~ 4 but with an extra decrease for early type galaxies which tend to 
have Xf\ increasing with A) while the dimensionless surface density k oc 1/S crit is a strongly 
increasing function of lens redshift for z\ <C z s . These trends are quantified in figure [13]. For early 
type galaxies, and in the /-band, the ratio of £;/k falls by about a factor of 3 as the lens redshift is 
varied from 0.2 to 0.4. Note also that high redshift clusters, if they exist, could give quite a strong 
lensing effect while being relatively inconspicuous in photometric surveys (the k curve is rather 
flat, while the Si >0 b s curve plunges rapidly). For example, a cluster at redshift unity produces a 
shear per surface brightness about one order of magnitude larger than the same cluster at redshift 
0.42 (this is neglecting the evolution of the cluster galaxies which tends to act in the opposite 
direction). 

To make an accurate prediction for the dimensionless surface density k then, it is not 
sufficient just to use the net surface brightness; one needs some idea of the distances to the 
galaxies. Ideally, one would obtain redshifts for all of the bright galaxies and then one could 
correctly account for this, but as yet we do not have this information. A solution to this problem 
is suggested by figure [5|. If we only use galaxies redder than V — I ~ 1.6 then we should 
see elliptical galaxies at z ^ 0.2, and later type galaxies only at z ^ 0.4. For an elliptical 
galaxy, we can read off the redshift from the color using figure |5] (we model the color-redshift 
relation as a linear trend fixed by the values at z = 0.2 and z = 0.4). We can then compute 
the physical luminosity, and thereby the mass (assuming a constant M/L), and thereby compute 
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SURFACE BRIGHTNESS AND SURFACE DENSITY vs REDSHIFT 
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Fig. 14. — The filled symbols show the observer-frame surface brightness (in the UH8K /-band) 
of a cluster lens as a function of lens redshift and for various galaxy types, assuming no evolution 
in the rest-frame surface brightness. Symbols encode the galaxy type as follows: EO = stars; Scd 
= squares; Sbc = triangles; Im = circles. The lines show the inverse critical surface density, which 
for a lens of a given physical surface density density, is proportional to the observed dimensionless 
surface density k, for sources at redshifts z s = (1.5, 2.0, 3.0). 
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it's contribution to the k field. For a single galaxy with observed magnitude m, and for an EdS 
cosmology, the contribution to an image of k from a galaxy which falls in a pixel with solid angle 
dfl is 

KdSl = - 4 ! G!Mq " (1 + ^)3^ilQO.4(M B0 -^)) (n) 

L <r(10pc) 2 uj s 



where M/ L is in solar units and to — 1 = 1/ vl + ^- 

If we generate a surface density prediction k in this manner we will, to a very good 
approximation have isolated the contribution from the early type galaxies at z ^ 0.2 alone. To 
see why, consider first the blue galaxies associated with the clusters at z ~ 0.42, which we have 
seen contribute at most about 30% of the cluster 5-band light. Some fraction of these — the very 
reddest — will survive the color limit of V — I — 1.6, but by the above prescription they will be 
assigned a redshift of ~ 0.2, rather than their actual redshift of 0.42 and so their contribution to 
k as calculated above will be suppressed in proportion to the lens redshift dependent factors in 
equation (fTT|): (1 + z) 3 ui(zi)ui s (zi, z s )10~ OAk ^ or about a factor 6 in this example, so, even if 
we assume that say 30% survive the color selection, the contribution to k will only be about 1/20 
of the true contribution of all the blue galaxies, or about 1/60 of the contribution of the elliptical 
galaxies. At higher redshift a larger fraction of later type galaxies survive the color selection, 
but for these the contribution to k will be suppressed by an even larger factor. Moreover, these 
galaxies are found to be very nearly uniformly distributed on the sky; thus their contribution to 
the fluctuating part of k — which is entirely responsible for the shear we are trying to predict — 
is expected to be very small and it should be safe to neglect them as compared to the early type 
galaxies at lower redshift. Also, our linear extrapolation of the color-redshift relation exceeds the 
actual predicted colors for z ^ 0.6, and this effectively suppresses the contribution from higher 
redshifts. 

With this color-based redshift estimation we obtain a prediction for k which will be correct 
if early type galaxies trace the mass with some constant M/L, (though missing some or all 
of the contribution from very distant structures where our selection function cuts off around 
or brighter). We should stress that this is a rather unconventional hypothesis. The general 
picture of the structure of superclusters, supported by the pressler (1980| ) morphology density 



relation, galaxy correlation functions ( |Loveday et al. 19^5| ) and studies of individual structures 



like the Virgo Supercluster, is of the elliptical galaxies residing in dense compact clusters, with 
the later type galaxies (which dominate the total luminosity) being much more extended. If this 
supercluster conforms to the norm, then the mass should be more extended than the early type 
galaxies and we should see this in the gravitational shear field, and if one incorporates statistical 
'biasing' (e.g. |Bardeen et al. 1986| ) then the mass would be expected to be still more extended. 



The resulting predicted k field smoothed with a Gaussian kernel of scale length 31" is shown 
in figure O. We show both the k just due to the z ~ 0.42 cluster galaxies and also for a broader 
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color selection. These were made with Mj L = 200/i which, as we shall see, is similar to the 
value preferred by our shear measurements. The differences are not large, though there is some 
contribution from the structure extending West from CL-S which lie at z ~ 0.3, and there are 
several other weaker features. 

The luminosity density field we measure is obtained from fluxes measured within apertures 
around the galaxies (these were taken to be three times the Gaussian scale length of the detection 
filter) and so will differ slightly from the total luminosity dependence because of two competing 
effects; the apertures are finite and so miss some of the light, but also sometimes overlap so that 
light gets double counted. To test how large the net bias is we have made estimates of the total 
extra-galactic surface brightness field in two ways: by summing the fluxes from the catalog and 
by using all the light in an image from which we have masked out all the stars. The results agree 
to better than 10 percent, so we conclude that we are measuring essentially the total luminosity. 

PREDICTED t(r) FROM CLUSTER E— GALAXIES PREDICTED ir(r) FROM ALL E— GALAXIES 
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Fig. 15. — Predicted dimensionless surface density k. Images have been smoothed with a 31" 
arcmin Gaussian low pass filter. The mean has been subtracted from each image. The left panel 
shows the predicted n only from galaxies in a narrow band of color around V — I ~ 2.2 for 
elliptical galaxies z = 0.42. The right hand panel shows the results for all galaxies redder than 
V — I — 1.6 and with contribution given by equation [TTJ. 



4.3. Predicted and Observed Shear 

Armed with the predicted k{t) from the previous section we can compute the shear field 
and compare with what we observe. To generate the predicted shear from the k images we solve 
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the 2-dimensional Poisson's equation by FFT techniques (padding the predicted k image out to 
twice its original size with zeros to suppress the effect of periodic boundary conditions) to obtain 
the surface potential $, and then differentiate this to compute the shear j a = M a i m $>i m /2. This 
provides us with an image of the shear field, and with this we can generate for each galaxy in our 
faint galaxy catalog the value of the predicted shear, and we can then ask, for instance, what is 
the multiple of the predicted shear which, when subtracted from the observed shear minimizes 
the summed squared residuals. We can also bin the galaxies according to predicted shear — here 
we treat the two components of the shear for each galaxy as independent measurements — and 
plot the average measured shear versus that predicted. To minimize the effect of shear generated 
by structures outside the field in this analysis, we subtract, from both the observed and predicted 
shears, the mean shear. 

If we assume that light traces mass and that the errors are predominantly in 
the shear measurement then to obtain the mass-to-light ratio we should compute 
a = Y,l° bs " 7 pred / J]7 pred ■ ^/P red 5 and then multiply this by our nominal value to obtain 
M/L B = 200/ia. We find a = 1.20, with this value giving a reduction in x 2 = XX 7 — Tmodd) 2 / ^ 
of 79 relative to the value for a = 0, from which we infer M/L B ~ 240/i. The correlation 
between the predicted and observed shear is shown in figure [16] and is significant at about the 
~ 10-sigma level. 

A weakness of this method is that the light was smoothed, whereas the shear values are 
unsmoothed. If light really traces mass fairly on all scales, then this would cause our estimate 
to underestimate the true mass to light ratio by a factor / d 2 kP(k)cu(k) 2 / J d 2 kP(k)u(k) where 
uj(k) is the smoothing filter and P(u) is the luminosity power spectrum. If mass were distributed 
around galaxies in halos of shape proportional to co(r), however, then we would obtain the correct 
Mj L. Conversely, if it were the case that the dark matter distribution is more extended than our 
smoothing filter then we will have underestimated the global M/L. We will return to this issue 
shortly. 

4.4. 2-D Surface Density Reconstruction 

We now show 2-D mass-reconstructions and compare with the projected luminosity density. 
We will also compare the results obtained using the V- and /-band observations separately, and 
with the level of noise expected due to the random intrinsic galaxy shapes. We use the Kaiser & 
Squires (1993) reconstruction algorithm. This is a stable and fast reconstruction method which 
has very simply defined noise properties; essentially Gaussian white noise. Its main drawbacks 
are that it does not properly account for patchy data and the finite data boundaries. We have found 
from simulations that the former is not a serious problem for these data, but the latter is worrying, 
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Fig. 16. — Predicted vs observed shear assuming M/L B = 200/i. Error estimates were obtained 
from an ensemble of catalogs with randomized shear values. 
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particularly for cluster CL-E, which lies right at the edge of the field and so some of the shear 
signal is lost, resulting in a suppression of the recovered mass. To get around this we proceed as 
follows: first we make a shear field image prediction from the constant M/L k prediction. We 
then sample this at the actual positions of our faint galaxies to generate a synthetic catalog (that 
which would have been observed with no intrinsic random shape or measurement noise), and 
then generate a reconstruction from that synthetic catalog which will have the same finite-field 
bias as the real reconstructions. To match the spatial resolution to that of the real reconstructions 
(a 45" Gaussian) we make the predicted shear with a smoothing scale 45"/v / 2 and make the 
reconstruction from the synthetic catalog with the same smoothing. Note that while correctly 
accounting for the finite field effect on structures within the field, the actual shear may still feel 
some effect from structures outside of the field. 

In figure [FT] we show the mass reconstruction from our combined faint galaxy sample on the 
left and that predicted (for a nominal M/ L B = 200h) on the right. There is generally quite good 
agreement between the mass and the light. All three of the main z ~ 0.42 clusters are clearly 
detected in the reconstruction, and the locations of these peaks coincide very accurately. We see 
no obvious indication that the mass in these clusters is more extended than the elliptical galaxy 
light. The lower redshift structure extending to the South-West of CL-S is also recovered in our 
reconstructions. There is some indication of a 'bridge' linking CL-S, CL-N, which is suggestive 
of the kind of features predicted by the 'cosmic-web' theory ( Pond, Kofman, & Pogosyan 1996| ). 
Part of this bridge is a clump at a ~ 3 /l 05 m 15", 5 ~ 17°16' which roughly matches a predicted, 
though somewhat weaker, feature. However, this is not at the supercluster redshift, and so part of 
the 'bridge' may be a coincidental projection. There is a rather conspicuous mass clump around 
a ~ 3 h 05 m 57", 5 ~ 17° 14' which has no apparent counterpart in the predicted field. 

There is considerable noise in the reconstructions. The noise power spectrum is white, while 
the galaxy power spectrum is relatively red, so we might expect to see better agreement if we 
further smooth both the light and mass maps, and to this end we show in the lower pair of panels 
of figure |T7| the results for a smoothing scale of ~ 90". 

Figure [18] shows mass reconstructions made from the I-band and V-band observations 
separately. These show generally very good agreement, with most of the features described above 
being visible in both images. The anomalous 'dark clump' lying between CL-S and CL-E is 
however not very evident in the V^-band reconstruction. Also shown is a reconstruction obtained 
using a catalog with the real faint galaxy positions but with the measured shear values shuffled 
(i.e. reassigned to the galaxies in random order). The fluctuations in this random reconstruction 
show the expected noise level due to random intrinsic ellipticities. We have made an ensemble of 
32 such realizations and we use these below to quantify the significance of our results. 
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Fig. 17. — Left panel shows the recovered mass using the |Kaiser & Squires (1993D reconstruction 
method. Right panel shows the predicted 2-D reconstruction using the shear predicted from the 
light-map (but sampled at the actual locations of the faint galaxies) which should have the same 
finite field biases as the real reconstructions. The lower panels show the same images smoothed 
with a Gaussian of scale 90" 
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Fig. 18. — Upper panels show the reconstructions from catalogs made from the I and V band 
images separately and show generally very good agreement. Lower left panel shows the early- 
galaxy M/L = 200/i prediction. Lower right panel shows the reconstruction from a catalog with 
randomized ellipticities, which indicate the expected noise fluctuations due to intrinsic random 
galaxy shapes. 
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4.5. Cluster Masses and Mass to Light Ratio's 

We now estimate the masses of the three z ~ 0.42 clusters. We will use these to derive 
mass-to-light ratios, and also to compare with independent mass estimates from the X-ray 
temperature and velocity dispersions. We compute the masses within circular apertures of radius 
1'.5 and 3'.0 from the reconstructions (which have been smoothed with a 45" Gaussian), and 
we estimate mass-to-light ratios by dividing the mass by the predicted mass within the same 
apertures. Thus, the M/L's we obtain should be corrected for finite field effects. The results are 
tabulated in table |3|. 



cluster 


r ap /arcmin 


K 


V 


h- l M/L 


hM/10 12 M Q 


CL-S 


1.5 


0.0318 ±0.0096 


3.3 


241 ± 72 


37.9 ± 11.4 


CL-E 


1.5 


0.0455 ±0.0164 


2.8 


182 ±65 


54.2 ± 19.6 


CL-N 


1.5 


0.0758 ±0.0100 


7.6 


342 ± 44 


90.5 ±11.9 


E+N+S 


1.5 


0.1586 ±0.0234 


6.8 


263 ± 38 


189.1 ±27.9 


CL-S 


3.0 


0.0131 ±0.0066 


2.0 


296 ± 149 


62.4 ±31.5 


CL-E 


3.0 


0.0206 ±0.0065 


3.1 


193 ±61 


98.1 ±31.2 


CL-N 


3.0 


0.0434 ±0.0081 


5.3 


466 ± 87 


206.9 ±38.7 


E+N+S 


3.0 


0.0764 ±0.0140 


5.4 


314 ±57 


364.7 ±66.9 



Table 3: Cluster masses and mass to light ratios. The mean dimensionless surface density K, the 
mass M, and the mass to light ratio have been corrected for finite field effects. The v value is the 
significance of the measurement determined by comparing with the rms mass measured for the 
ensemble of randomised catalogs. 



Table |3] shows an increase in Mj L with aperture radius for all three clusters. For CL-E, CL-S, 
this is slight and not particularly significant. For CL-N the increase is larger, but we caution that 
there is a foreground cluster which lies just to the East of CL-N and which falls within the larger 
aperture. While we have argued above that there should in general be a rather weak contribution 
from low redshift structures, in this case the foreground cluster is particularly bright, and so 
may contribute significantly to the increase in Mj L for CL-N. Redshift information is needed 
to properly quantify this. We note that these M/L values, averaging around 260/i for our 1'.5 
apertures (which give the highest signal to noise for Mj L) are much lower than has been found 
for other, generally more massive, clusters like A 1689 ( |Fischer & Tyson 1997| ). In so far as only 
a tiny fraction of the Universe resides in such super-rich systems, it may well be that the value 
obtained here is more representative. It is also interesting to note that these Mj L's only include 
the contribution from early type galaxy light in a rather narrow band of color. CL-N was found 
by pabricant, Bautz, & McClintock (19941 ) to contain a number of bluer 'Balmer-line' galaxies, 
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about ~ 25% of the total, which would bring the total M/L is more in line with that for CL-S, 
which they found to be almost entirely early galaxy dominated. 

It is of interest to compare these masses with those obtained from the CL-S temperature and 
with velocity dispersions for all 3 clusters. For CL-S, our mass is in good agreement with the 
X-ray mass (which in turn is in good agreement with the velocity dispersion). For CL-N, we 
model the galaxies as an isothermal gas with profile index dlogn/dlogr ~ —1.5 as measured 
from the observed counts profile to find a predicted k ~ 0.93 (for the 1'.5 aperture) as compared 
to the measured value 7t = 0.0 76 ±0.01, again in reasonably good agreement given the statistical 
and systematic modeling uncertainties. For CL-E the same exercise gives a predicted surface 
density about a factor two higher than what we measure, but this may be due to departures from 
assumed equilibrium in this clearly non-relaxed system. 



4.6. Light-Mass Cross- Correlation 

In this section we show the results of cross-correlating the mass and the light. Our goals are 
both to determine the Mj L parameter and also to test the hypothesis of a constant Mj L which is 
independent of scale. We perform the correlation both in real space and Fourier space. 



The real space cross -correlation function is shown in figure [I9|. In computing this we 
padded the source images with zeros to twice the original image size. Also shown is a typical 
random realization from the ensemble used to estimate the noise fluctuations. We see a strong 
cross-correlation peak at close to zero lag. The significance, estimated as the strength of the 
zero-lag correlation relative to the rms found from our ensemble of randomized catalogs is 
9.3-sigma. The correlation strength at zero lag implies a mass-to-light ratio M/L = (224 ± 24)/t 
in solar units. 

To see if Mj L changes with scale, in figure |2(] we show the profile of the luminosity-mass 
cross-correlation and the luminosity auto-correlation function. These are very similar in shape. 
We caution that these should in no way be taken to be indicative of the mass-luminosity 
auto-correlation function in general; the field is small, so the 'cosmic-variance' or sampling 
variance is very large, and the field is also unusual in that it was chosen because it contains a 
prominent supercluster. What this figure shows though is that aside from a slight enhancement 
of the luminosity auto-correlation function at very small lag — which might plausibly be due 
to dynamical segregation — the cross- and auto-correlation functions have very similar shapes 
indeed. This suggests that on scales ^ 1' the early type galaxies trace the mass quite faithfully. 

Another way to see if M / L changes with scale is to perform the correlation in Fourier space. 
The discrete Fourier transform of the k image provides a set of estimates of the actual mass 
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Fig. 19. — Cross-correlation of light with mass reconstruction and with randomized catalog 
reconstruction. Left panel shows that actual cross correlation and right panel shows the result 
of cross correlating the light with a k reconstruction generated from a catalog with randomized 
(shuffled) shear values. 

transform which have essentially statistically independent measurement errors. The errors are 
exactly independent only if the data fully cover the field - here there are some gaps in coverage 
like the N-W corner which introduce slight correlations between neighboring Fourier coefficients. 
This slight covariance is properly accounted for by our Monte-Carlo error estimation. The results 
are shown in figure [H| Here the transform was taken without zero padding. If we use all modes 
with k < 8 times the fundamental frequency for our square box (corresponding to wavelengths 
A = 0.75 — 6/i -1 Mpc) we obtain a correlation which is significant at the the 8.3-sigma level and 
we obtain Mj hL = 233 ± 28. If we split the modes into low and high frequency subsets we find 
M/hL = 280 ± 42 for < k < 4 and M/hL = 176 ± 38 for 4 < k < 8 so again, we find a slight, 
though barely significant, increase in Mj L with scale. 



5. Discussion 

The foregoing analysis shows in a number of ways that we have clearly detected the 
gravitational shear from mass concentrations in this field. All of the major concentrations apparent 
in the X-ray and optical images are detected in the mass reconstructions, as well as a foreground 
structure at z ~ 0.3. While the mass reconstructions are somewhat noisy, we have shown by a 
number of different (though not independent) methods that there is a strong correlation between 
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Fig. 20. — Azimuthally averaged profile of the mass-luminosity cross -correlation function (boxes 
with error bars) and the luminosity auto-correlation function (circles). In normalizing these we 
have adopted a mass-to light ratio of Mj L = 280/i, somewhat larger than that obtained by simply 
comparing the correlation strengths at zero lag. 
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Fig. 21. — Cross-correlation of light with mass reconstruction in Fourier space. Left panel 
shows that actual cross correlation and right panel shows the result of cross correlating the light 
with a k reconstruction generated from a catalog with randomized (shuffled) shear values. Box 
symbols are for modes with A = 1.5 — 6/i _1 Mpc and + symbols are for modes with wavelength 
A = 0.75 - l^/i^Mpc. 
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the observed shear and that predicted assuming early type galaxy light traces the mass, which 
allow fairly precise estimate of the total mass associated with these galaxies, and we find that the 
Mj L parameter does not vary much with scale. 

The dynamical state of the z ~ 0.42 super-cluster is that most of the mass in this region is 
concentrated into the three clusters. Following [Fabricant, Bautz, & McClintock (1994D we have 
modeled CL-S as a test particle on a radial orbit around CL-N and CL-E as a test particle on a 
radial orbit under the combined attraction of CL-N and CL-S. We find generally similar results, 
but as our masses our somewhat lower than they assumed (by scaling to Coma using the velocity 
dispersion), our orbital solutions are less tightly bound. We find that CL-S is only marginally 
bound to CL-N. Equivalently, the mass of CL-N, if spread over a sphere or radius equal to the 
distance to CL-S, gives a density about equal to that of a critical density Universe at that epoch. 
We find that CL-E is on an unbound trajectory, though the latter conclusion could possibly be 
modified if it turns out that CL-N, CL-S have massive neighbors outside of the field we have 
surveyed. While the clusters are bound and have collapsed, the super-cluster as a whole is still 
rapidly expanding and may never turn around. 

As remarked earlier, the conclusion that early type galaxies trace the mass is rather surprising. 
The conventional picture of the morphological segregation in superclusters is heavily influenced 
by the Local Supercluster. There one finds flBinggeli, Tammann, & Sandage 19871 ) a E-rich core of 
dimension ~ l/i _1 Mpc surrounded by a spiral rich 'halo' extending to large radii (~ 10/i _1 Mpc) 
roughly as 1/r 2 and which dominates the total light. This picture is supported in a statistical 
sense by the pressler (1980D morphology-density relation and also by galaxy correlation studies 



( Loveday et al. 199~5[ ). It is commonly assumed either that galaxies in general trace the mass, or, 



in biased theories and in hot dark matter models, that the mass distribution is even more extended. 
What we have found is in conflict with either of these pictures; the mass around the clusters in this 
field is no more extended than the early type galaxies. 

Exactly what this implies for the global density parameter Vt remains somewhat unclear 
because of the possibility of bias; as with other dynamical estimates of cluster masses one can 
estimate Mj L and then extrapolate to the entire Universe as a simple accounting exercise to 
obtain the mass density as the product of Mj L with the luminosity density measured from redshift 
surveys, but there is of course no guarantee that the cluster M/L's are really representative. 
Nonetheless, and with the foregoing caveat, it is of interest to make this extrapolation to obtain f2 
assuming that galaxy formation is unbiased and one can then say what level of bias is implied for 
other values of O. 

If we assume that galaxies of all types formed with their usual abundance in this region then 
we have to ask what has happened to the late types. One possibility is that they are still there and 
are more broadly distributed than the early types as is the case in other 'normal' systems (we 
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know from the spectroscopic studies that they are not in the clusters themselves). Our data do 
not really support this picture, but we cannot firmly rule it out either. If so, these galaxies have 
very little mass associated with them, and the total density of the Universe is essentially just the 
luminosity density in early type galaxies times our M / L and is very small. [Loveday et al. (199"2|) 



for example, find a total luminosity density of ~ 2 x 10 8 /iL Q /Mpc 3 of which about 20% is due 
to the early types, and we thereby find Q ~ 280 x 0.2/1500 ~ 0.04. The other possibility is that 
the late type galaxies in this region have, like the early types, fallen into the clusters, but have had 
their disks extinguished by ram pressure. If there were much mass associated with these galaxies 
then we should have seen very high mass-to-light ratios for these clusters, but we didn't, so again 
this would suggest there is very little mass associated with late types. However, the bulges of the 
spirals will survive, and would have inflated the cluster luminosity functions at the faint end and 
the true value of f2 is then slightly higher. 

If galaxy formation was unbiased then, our results suggest a very low value for the density 
parameter, and, if h is at the lower end of the allowed range, this is mostly or all baryonic. We 
do not know that early type galaxies are unbiased tracers, but the assumption that they are is at 
least consistent with the general idea that they formed early (so any statistical bias they may have 
had a the time of formation may have been washed out) and that they seem to be evolving slowly, 
and therefore seem to be obeying the continuity equation in an average sense. To avoid this rather 
radical conclusion one would have to invoke a bias. However, it would seem to be very hard to 
accommodate these data with say a flat EdS model, since in that case one has a region with net 
mass density at about or even slightly below the mean, but where the formation of early type 
galaxies was apparently biased upwards by a factor 20 or so. This seems unreasonable. More 
modest positive bias cannot be ruled out, and it is also possible that the early type galaxies are 
'anti -biased'; i.e. the abundance of these galaxies per unit mass is lower than average in clusters, 
in which case the global fl is still lower than the above estimate. 

We note that our mass determination is blind to any uniform density component, so we 
cannot rule out the possibility of a contribution to £1 from very hot dark matter or an ultra-light 
scalar field, the latter having an effective Jeans length on the order of the geometric mean of the 
horizon size and the Compton wavelength and which could be very large. 

While the density parameter we find seems low as compared to other estimates from 
e.g. virial analysis ( Pavis & Peebles 1983| ; |Carlberg et al. 199^| ), this is not because our cluster 



M/L values are much lower, it is because we do not extrapolate to the Universe at large assigning 
the same Mj L to late type galaxies. The interesting new conclusion from our results is that if late 
types do have the same net M / L as early types — sofi ~ 0.2 — 0.3 — then their formation in 
the rather large region of space originally occupied by the matter now in the supercluster must 
have been strongly suppressed and the galaxy formation process must have arranged to create 
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essentially only early type galaxies, and with an abundance several times higher than average to 
compensate. 

Our cluster MJ Ls are considerably lower than that obtained from weak lensing analyses 
of super-massive clusters like A1689 (e.g. |Fischer & Tyson 1997| ), and it is interesting to note 
that the most massive of the three clusters here has a higher M/L than the other two. If these 
results are indications of a general trend then this would seem to argue against a positive bias for 
clusters, since the most extreme clusters seem to be less positively biased than more common 
poorer clusters like those studied here. Equally if f2 is really 0.2 — 0.3 as often assumed then 
this requires non-monotonicity of bias with both extremely massive clusters and field galaxies 
being essentially unbiased, but with intermediate mass clusters being positively biased. This 
seems contrived. Finally, it is also interesting to briefly compare our results with peculiar velocity 
analyses which provide an alternative probe of the total (dark plus light) mass distribution. Barring 
some very hot or ultra-light mass component, our low would seem to be very hard to reconcile 



with the generally high values found from large-scale 'bulk-flows' ( pekel 1994| ). On the other 
hand, if there is really very little mass outside of clusters, as our data suggest, then this would be 
compatible with the very cold nature of the Hubble flow in the field (|Sandage 1986|; Brown & 



Peebles 1987; Pstriker & Suto 1990^ 7 
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